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Abstract 

We consider the low energy limit of a stack of N M-branes at finite temperature. In this limit, 
the M-branes are well described, via the AdS/CFT correspondence, in terms of classical so- 
lutions to the eleven dimensional supergravity equations of motion. We calculate Minkowski 
space two-point functions on these M-branes in the long-distance, low-frequency limit, i.e. 
the hydrodynamic limit, using the prescription of Son and Starinets |hep-th/0205051 . From 
these Green's functions for the R-currents and for components of the stress-energy tensor, we 
extract two kinds of diffusion constant and a viscosity. The N dependence of these physical 
quantities may help lead to a better understanding of M-branes. 
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1 Introduction 



The interacting, superconformal field theories (SCFT) living on a stack of N M2- or M5- 
branes are not well understood. An improved understanding of these M-branes should lead 
eventually to a better understanding of M-theory itself, a theory that encompasses all the 
different super string theories and is one of the best hopes for a quantum theory of gravity. 
While the full M-brane theories remain mysterious, the low energy, large N behavior is 
conjectured to be described well, via the AdS/CFT correspondence PQ 12 El, by certain 
classical solutions to eleven dimensional supergravity equations of motion. Recent work on 
AdS/CFT correspondence by Son, Starinets, and Policastro (HE] provides a prescription for 
writing Minkowski space two-point functions for these types of theories. We take advantage 
of this prescription to calculate viscosities and diffusion constants for M-brane theories in 
this low energy limit, thus generalizing the work of |Hj for D3-branes. 

Policastro, Son, and Starinets |5 3 used their Minkowski space prescription to investigate 
the low-frequency, long-distance, finite temperature regime of the D3-brane theory. There is 
lore that the long-distance, low-frequency behavior of any interacting theory at finite temper- 
ature can be described well by fluid mechanics (hydrodynamics) |0j. Although not rigorously 
proven, the idea is well supported by physical intuition about macroscopic systems. Hydro- 
dynamics in turn provides rigorous constraints for the form of Minkowski space correlation 
functions. Once a few viscosities, diffusion constants, and other transport coefficients are 
known, the two-point functions are completely fixed [Zj. 

Indeed, the authors of [5] found that the form of the Green's functions calculated from 
supergravity was completely consistent with hydrodynamics for this D3-brane theory. More- 
over, from these Green's functions, they were able to extract a transport coefficient (a shear 
viscosity) and diffusion constants. The authors' prescription for Minkowski space Green's 
functions is a modification of the prescription for calculating Euclidean space Green's func- 
tion developed in [2JE1- Both prescriptions use the gravitational description of the low energy 
theory on the brane for calculating correlation functions. 

It is true that these same Green's functions on a D3-brane can be calculated in a more 
traditional way using weakly coupled gauge theory |Hj. The low energy theory living on a 
stack of iV D3-branes is described via the AdS/CFT correspondence, alternately as M = 4 
SU(N) super Yang Mills or as type IIB supergravity in a AdS§ x S 5 background. While 
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the gravity calculation is good at strong 't Hooft coupling A (large curvature), the gauge 
theory calculation is good only for small A. However, the N and temperature dependence 
of these transport coefficients and diffusion constants should be and indeed is universal and 
independent of the method of calculation. Remarkably, the gravity calculation appears to 
be technically simpler than the equivalent calculation for weakly coupled gauge theory. 

In contrast, the M-brane theories have no such alternate gauge theory description. How- 
ever, there are still low-energy supergravity solutions from which we can extract analogous 
Minkowski space Green's functions. Moreover, we can hope that M-branes, like D3-branes, 
have a hydrodynamic regime. 



Indeed, based on the Minkowski space prescription of Son and Starinets £fj . in the low- 
frequency, long-distance, finite temperature regime, we find that these M-brane theories have 
two-point functions which are completely consistent with a hydrodynamic interpretation. 
Generalizing jH] , we calculate two-point functions for the conserved R-symmetry current and 
for components of the stress-energy tensor for these M-branes. Moreover, from these Green's 
functions we are able to extract corresponding diffusion constants and also the viscosity. The 
N dependence of these physical quantities may lead to a better understanding of M-brane 
theories. 

It should be emphasized that the Son and Starinets prescription 4J for calculating 
Minkowski space correlators is not completely justified, and that this paper provides some 
limited additional evidence for their prescription. Turning the logic of the last two para- 
graphs around, we can argue that on general grounds we expected the M-branes to have 
a hydrodynamic description. Thus, it is reassuring that the prescription of [I] does indeed 
produce Green's functions with the appropriate behavior. Also, we get results which are 
internally consistent. On general hydrodynamic grounds, we expect that D = 77/(6 + P) 
where D is the diffusion constant calculated from the stress energy tensor, 77 the viscosity, 
e the energy density, and P the pressure. We check that this equation holds both for the 
M2-branes and M5-branes. 

We begin by reviewing some essential facts about the non-extremal M2- and M5-brane 
backgrounds. As we are working at finite temperature, the extremal AdS$ x S 7 and AdS-j x S 4 
supergravity solutions are not adequate. We need their nonextremal generalizations where 
we can associate a Hawking temperature to the horizon. 
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Next, we consider R-current correlation functions. Both M2- and M5-brane supergravity 
solutions have an R-symmetry which one can think of roughly as the rotational symmetry 
of the transverse sphere. From the form of the thermal R-current two-point functions, we 
extract a corresponding R-charge diffusion constant. 

Having warmed up with the R-current correlators, we proceed to the more complicated 
example of stress-energy tensor two-point functions. From components of these two-point 
functions, we extract a diffusion constant and a viscosity for each M-brane theory. 

Finally, we end with some comments about the iV dependence or lack thereof of the 
various diffusion constants and transport coefficients calculated. The motivation for this 
work came in large part out of the hope that some of these N dependences might shed light 
on the underlying M-brane theories. 

2 The Nonextremal Eleven Dimensional Supergravity 
Backgrounds 

We begin by reviewing some essential facts about the non-extremal M2- and M5-brane 
supergravity solutions. Roughly speaking, these solutions represent what happens when a 
stack of M-branes is placed in flat 11-dimensional space and given some finite temperature 
T. The 11-dimensional space, close to the M-branes, separates into a product of a sphere 
and an asymptotically anti-de Sitter space. A horizon with Hawking temperature T forms. 
These supergravity solutions are solutions to the equations of motion following from the 
eleven dimensional supergravity action [5] 

-V / d n x(-g) 1/2 K _ -L / (> 4 A *F 4 + \a z A F 4 A F 4 
2kJi J 4:K 2 n J \ 3 

where «n is the gravitational coupling strength, and dA 3 = F 4 . 
2.1 M5-brane 

For the M5-brane, the nonextremal metric is 

ds 2 = H(r)~ 1/3 [-f(r)dt 2 + dx 2 ] + H(r) 2/3 



dr2 +r 2 dnj 



where H(r) = l + i? 3 /r 3 and f(r) = 1— r^/r 3 . The quantity dx 2 is a metric on flat, Euclidean 
IR 5 . The term dQ 2 is the metric on a unit four sphere S 4 . The four form flux F 4 from the 
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M5-branes threads this S 4 : 

F 4 = 3i? 3 vol(S 4 ). (3) 

The quantization condition on the flux implies that N 3 ^ = 2 7 n 5 R 9 where N is the number 
of M5-branes and Kn is the eleven dimensional gravitational coupling strength [TU] . 

Taking the near- horizon limit r -C R, according to the AdS / CFT prescription QjJ El 
we can "zoom in" on the five-brane theory dynamics: 

ds 2 - ^ [-f{u)dt 2 + dx 2 ] + + R 2 dVt\ . (4) 

uR L J f(u) u z 

We have made the coordinate transformation u = r^jr. Spatial infinity, which is also now 
the boundary of an asymptotically anti-de Sitter space, corresponds to u — 0. There is a 
horizon at u = 1. The Hawking temperature of this horizon is 

3 r 1/2 

Another important quantity characterizing this supergravity solution is the entropy den- 
sity, which one finds by multiplying the horizon area by 27r/Acf 1 and dividing out by the 
volume of the gauge theory directions, x 1 , x 2 , . . . , x 5 (TTJ: 

2 7 7T 3 

S = -y-N 3 T 5 . (6) 

We will use index conventions where fi,u, . . . refer to the asymptotically AdS directions, 
a,{3, . . . index the M-brane directions only, and i,j, . . . index the spatial M-brane directions. 

2.2 M2-brane 

An analogous nonextremal supergravity solution exists for a stack of M2-branes in eleven 
dimensional space. Now the metric takes the form 

dr 2 

W: 

where H(r) = 1 + i? 6 /r 6 and f(r) = 1 — r^/r 6 . The four form field strength is easier to think 
of in a dual language as a seven form field strength: 

*F 4 = F 7 = QR\o\(S 7 ) . (8) 

The quantization condition ^U] on the field strength reveals that i? 9 7r 5 = N 3 / 2 k\ 1 \/2. 



ds 2 = H{r)- 2 ' 3 [~f{r)dt 2 + dx 2 + dy 2 ] + if(r) 1 / 3 



+ r 2 dn 2 



(7) 
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We can again take a near horizon limit, r -C R, to find 

ds 2 - 4m h/(^ 2 + ^ + d V 2 ] + TT^T^T + ^ > ( 9 ) 
where w = r 2 /r 2 . The Hawking temperature is 

T = — ^ (10) 

27T.R 3 ' 1 J 

and the entropy density is [11 J 

27 v ; 

3 R-charge Diffusion for M2- and M5-branes 

The R-charge interactions are mediated in the bulk by a gauge field F*. The starting point 
for calculating two-point R-charge correlators is the usual Maxwell action in the nonextremal 
backgrounds given above: 

1 



S = -jj- / d d x^F^ a . (12) 

Here d is equal to 4 for the M2-branes and 7 for the M5-branes. The tensor is a metric on 
the asymptotically anti-de Sitter space. The calculation is similar to the Euclidean calcula- 
tions in ^2]; however we work with a Lorentzian signature and in a nonextremal background. 
The constant gsc can be set by compactifying the eleven dimensional supergravity action 
(JTJ) on a S 4 or S 7 . We will ignore this overall normalization for the Green's functions for 
now. The reason is that our main interest in this section is the diffusion coefficient for the 
R-current which can be obtained simply from the location of the pole in the corresponding 
retarded Green's functions. The overall normalization of the Green's function is also related 
to the diffusion constant, but in a more complicated way via a Kubo-type formula. 

Our calculations closely follow [5 1. We work in the gauge A u = 0. We use a Fourier 



\= I ^k^ M ^M^)- (is) 



decomposition 

A = ( 

(270^ 

Rotational invariance in the spatial directions allows one to simplify things further by choos- 
ing q° = u, q 1 = q, and all other q l = 0. The equations of motion for the are 

1 --9u W=~gg w g v ° {d P A a - d a A p )] = o . (14) 



At this point, it becomes convenient to analyze the M2- and M5-brane cases separately. 
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3.1 R-charge and M5-branes 



The equations of motion (|T3j) for reduce to 



+ fq 5 A' x 

K ~-K-—f {^bQbAc + q\A t ) 
u uj v 

A" x - -A' x + j -A' x + \ (q 5 u 5 A t + u 2 A x ) 
u j uj z 

"> J "J 



where t is x°, x is x 1 , and a stands for any of the other x l . We have introduced 



, (15) 
, (16) 

, (17) 

, (18) 

" 5 =4^ ' q ^^T q - (19) 
The prime denotes the derivative with respect to u. The peculiar combinations of g|, ^5^5, 
and q 2 in the equations guarantee gauge invariance under the residual transformation A t — ► 
At — cuA and — > A x + qA. The first three equations are dependent; equations (fT5j) and 
(fT5|l imply equation (JTTj). 

Because of our choice of g-vector, the R-charge diffusion appears only in the A t and A x 
sector, and we begin with the first two equations, (fT5|) and (|16)). These two equations can 
be combined to yield a single equation for A' t : 



f_ A" 

f 



A'!' + J -A'l + 



vlf ~ ff A , 



up 



A =0 . 



(20) 



This second order equation for A' t does not appear to be analytically tractable. However, a 
solution can be obtained perturbatively for small q^ and tv^. Following we determine the 
behavior of A' t near the singular point u — 1. Substituting A[ = (1 — u) a F(u), where F(u) 
is a regular function, one finds that a 2 = —uj 2 /9. The "incoming wave" boundary condition 
described in |E] forces us to choose a = —iu 5 /3. 

Next, we solve for F(u) perturbatively: 

F{u) = F (u) + L0 5 Ft(u) + qlG x {u) + 0{ul L0 5 q 2 5 , g 5 4 ) , 



where F (u) = uC and 



F x {u) 
GAu) 



iC 
C 



u 



[u 



l) + g/l(«) 



u 

71 



1 / s U „ , s 

-(1 -«) + _/„(«) 



(21) 

(22) 
(23) 
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where 

h(u) = In ; f 2 (u) = tan [—^~ J ~ 3 ■ ( 24 ) 

One of the two integration constants for Fq(u) is set by requiring that Fq(u) is well-behaved 
at u = 1. The integration constants for all higher F{(u) and Gi(u) are set by requiring that 
these functions vanish at u = 1 as well. 

The constant C can be related to the boundary values A® and A x using (fTBJ): 

C= ^ A °\"l A < . (25) 
zo; 5 - 

The pole in C is the same pole that appears in the retarded Green's functions, as we will 
presently see. Having obtained A' t , A' x follows from lJI5jl . 
The solution to pSjl can be obtained in a similar fashion: 



A a = aK ' ' } + O(col co 5 ql 4) (26) 



where 



Ku) = 1 + f - - ti^fM ■ ( 27 ) 

The Green's functions can now be calculated from the terms in the action which contain 
two derivatives with respect to u: 

1 



(28) 



S = -—z- / dud b x V^gg^'dvAidvAj + ... 

•y SG 



fi 1 

, du d x — 

2&g 2 SG J u 



5 

/2 

x' 

i=l 



Recall from jl] the procedure for calculating these Minkowski space Green's functions for a 
scalar <f>(u). We extract the function A(u) that multiplies (d u <f)) 2 in the action: 

S = - [ dud & xA{u){d u (j))' 2 . (29) 



2 . 

Next, we express the bulk field 4> via its value 0o at the boundary u = 0, (f>(u, q) = f q (u)(po(q). 
By definition f q (0) = 1. Moreover, we impose an incoming-wave boundary condition on f q (u) 
at the horizon u = 1 (when g is timelike). The retarded (Minkowski space) Green's function 
is then defined to be 

G R (q) = A(u)U q (u)dJ q (u)\ u=0 . (30) 
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Using this prescription, we find that the retarded Green's functions for the R-current are 

G$ = C . q ° n , +■•• , (31) 
iuj - D R q 2 

G£ = G£ = -C " q8 * +... , (32) 
iuj - D R q 2 

oj 2 5 ab 

G a xx = C- — - H , (33) 

iuj - D R q 2 

Gt = -C6 ab (iu + 2D R q 2 ) + ■■■ . (34) 



where 

3/2 

C = (35) 

and 

Dr = -4^ • (36) 

As expected, the G aa Green's functions have no pole while the others do. From the location 
of the pole, we can read off the diffusion coefficient for the R-charge, D R . We regard this 
value of D R as a prediction for the theory living on a stack of M5-branes at finite temperature. 
The power of T in D R is forced by dimensional analysis. However, it is interesting that this 
value for D R is iV independent. 

This expression for D R is subject to two kinds of correction. First, the supergravity 
approximates these M-branes well only at large N. Thus, there could be 1/N corrections. 
The second correction is not really a correction to D R but to the location of the pole itself. 
We could easily calculate the to higher order in q 5 and u 5 . In this case, we would get 
corrections of order 0(q 2 ) to the location of the pole. 

Although we did not solve for g$G exactly, we can count powers of iV and T in C. The 
power of T is forced by dimensional analysis. Tracing powers of iV is relatively easy. The 
coupling g SG is essentially K\\ multiplied by lots of iV independent compactification factors, 
and /en ~ iV~ 3 / 2 . 

3.2 R-charge and M2-branes 

We now redo this same calculation in the nonextremal M2-brane background (J9jl . The 
differential equations for take the modified form 



u 2 A' t + q 2 fA' x = 0, (37) 



A" - {u 2 q 2 A x + q\ A t ) = 0, (38) 
K + L A ' x + ^L(uj 2 q 2 A t + u 2 2 A x ) = 0, (39) 



K + fK + Tn H-h!)A = 0. (40) 



— A' 

f y ' 4/ 

where x° = t, x 1 = x, and x 2 = y. We have defined the quantities 



" 2 "^t" ]q2 =2^T q - (41) 
This system is very similar to the one encountered in the previous section. Equations (|37j) 
and (jHEj) imply equation (jHSJl- We combine equations ((57)) and ((55)) to give a single differential 
equation for A J alone: 

A7 + ^ + ^(u£ - fq 2 2 )A[ = . (42) 

We make the substitution A' t = (1 — u)~ lul2 / 6 F(u) and solve for F(u) perturbatively in uj 2 
and q 2 : 

F{u) = C(l + w 2 Fi(u) + g 2 2 Gi( M ) + . . .) (43) 

where 

Fi(v) = ^h(u) + ^=f 2 (u), (44) 
Gi(u) = -l^hiu). (45) 

Using ((55)1. we can solve for C in terms of the boundary values of At and 

o; 2 g 2 Ag + q 2 2 A ° t 
2(iu 2 - \ql) 

The pole in C will be the same pole that appears in the R-current Green's functions and 
hence is related to the diffusion coefficient. 

The perturbative expression for A' x can be obtained from (J37j) . We have already done 
the necessary work for calculating A y . Note that (J4U)) is the same differential equation as 
flUJ). Thus A y = (1- u)~ luJ2/(i F(u) with F( u) given by (|43|). The only change is that now 
C = A y + 0(LU 2 ,q 2 ). 

To extract the Green's functions, we need to isolate the terms in the action with two u 
derivatives: 

S = ^ J dud 3 x [A? - fA' 2 - fA' y 2 + ...} . (47) 
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UJ 2 q 2 ST. x T q 2 sx t ( s 



From this expression, it is straightforward to see that 

g 2 6 ab 



^xt ~ ^tx 



G 



ab 

yy 



g 2 SG (iu - D R q 2 ) ' 
uq5 ab 

9sc( iuJ ~ D Rl 2 ) 
uo 2 5 ah 

g 2 SG (iuJ - D R q 2 ) ' 

5 ab 

=- (iu - D R q 2 ) 

9sg 



where 1/ggQ ~ N 3 ^ 2 and the diffusion coefficient is 

3 



D 



R 



AttT 



This expression for D R is again N independent. 



(48) 
(49) 
(50) 
(51) 



(52) 



4 Stress Energy Two Point Functions 

To obtain more diffusion constants for our M-brane theory in the hydrodynamic limit, we 
compute the two-point function of the stress-energy tensor. According to the AdS/CFT 
prescription, this two-point function is related to small perturbations of the metric in the 
bulk theory. In particular, we consider g^ v = gffl + h^ v where the unperturbed metric g^J 
is given by the asymptotically AdS pieces of the full eleven dimensional metrics Q and (@J. 
The transverse spherical parts of the eleven dimensional metrics we leave unchanged. 

For this calculation, we can ignore the spherical piece of the metric and the four form field 
strength F 4 . They do not enter at first order in the perturbations. For the asymptotically 
AdSj case, i.e. the five-brane case, we consider the compactified action 

i? 4 Vol(S 4 ) 



S 



2^ii 



J dud 6 x^g{TZ-2A) + 2 J d % x\[^g®K 



(53) 



In the above, 1Z is the Ricci scalar and A = — 15/4i? 2 is a cosmological constant arising from 
the integral over I-F4I 2 in the full eleven dimensional action. The metric g^g is the metric 
induced on the boundary u = while K is the extrinsic curvature of the boundary. There 
is an analogous action for the asymptotically AdS 4 case. 

R 7 Vo\(S 7 ) 



S 



du d 6 x^g(Jl - 2A) + 2 



(54) 
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where the cosmological constant is now A = —12/R 2 instead. 

The first step in computing the stress-energy two-point function is to solve the linearized 
Einstein's equations 

n $ = d^2 Ah ^ (55) 
where d is either 4 (for asymptotically AdS^) or 7 (for asymptotically AdS 7 ). By 1Z$, we 
mean everything in the Ricci curvature 1Z^ U that is linear in h^ v . 

The metric perturbations split up naturally into groups as follows. We take a Fourier 
decomposition of the metric perturbation, assuming that h^ v depends on t and x as e - luJt + l i x ^ 
defining x° = t, x 1 = x, x 2 = y, and so on. We choose a gauge such that A^ u = 0. For 
the M5-brane case, there is a rotation group SO {A) acting on the directions transverse to 
u, t, and x. The perturbations split according to whether h^ u is a tensor, a vector, or a 
scalar under these rotations. For the asymptotically AdS± case, there is only one direction, 
the y direction, transverse to the others. However, we can still consider the effect of sending 
y — > —y. Under this transformation, the metric perturbations split into two groups according 
to whether h^ v has an odd or an even number of y indices. Stretching definitions, we will 
call the metric perturbation with only one y index a vector mode. There is one remaining 
perturbation involving components of h^ v with both no and two y indices which we will call 
the scalar mode. 

The modes are useful in different ways. The tensor mode allows us to compute a shear 
viscosity of the boundary theory using a Kubo formula. There is a diffusion pole in the vector 
mode which will allow us to calculate a diffusion constant. There is also a Kubo formula 
for these vector modes which will allow for another check of the viscosity. The scalar mode 
describes sound propagation on the boundary. 

We will begin by computing the diffusion constant from the poles in the vector modes, 
both for asymptotically AdS± and AdS 7 . The calculation is very similar to that of the Re- 
charge diffusion constant. The relation should not be at all surprising because from an eleven 
dimensional standpoint, the A^ vector potential is is in part a metric perturbation of the 
form Hm^ where M is a spherical direction and /i is an asymptotically AdS direction. Next, 
we will analyze the tensor mode for asymptotically AdSj. We leave sound propagation for 
future work. 

Sound propagation is the most involved of the three types of metric fluctuations. To 
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consider a tensor fluctuation for the M5-brane, h yz is the only component of the metric 
fluctuations that needs to be nonzero. For the vector fluctuations, two components of 
need to be nonzero, for example h xy and ht y . For the scalar modes, all of the diagonal 
components of h a p plus h xt must be nonzero. The resulting system of differential equations 
is less tractable than for the tensor or vector modes. 

4.1 Vector Modes and a Diffusion Constant for M5-branes 

We consider a metric perturbation of the form h ty ^ and h xy 7^ with all the other h^ v = 0. 
Moreover, we make a Fourier decomposition such that 

hi = e - iuJt+iqx H x (u) , (56) 
h y = e -i"t+n* Ht (u) . (57) 

The linearized Einstein's equations for H x and H t are 

oj 5 H' t + q 5 H'J = 0, (58) 

H'^--H' t -—(u^H x + qtH t ) = 0, (59) 
u uf 



K-^—tK + ^-^Ht + uiiH,) = 0. (60) 
uf uf 2 



This system of equations is very similar to the systems (|T3j) - (fT7j) and ()37|) - (pl9~|) we solved in 
the previous sections and is tractable using exactly the same methods. We combine the first 
two equations (jB^j) and (jBTJj) to get a single equation for H' t : 

H? + 2J ^K + \M ~ f<& + ^f)H> = . (61) 
uf uf 2 

We make the ansatz H' t = C(l — u)~ tuJ5 / 3 F(u) and solve for the regular function F{u) 
perturbatively in co> 5 and q§: 

F(u) = u 2 + iu 5 Q(« 2 - 1) + l -u 2 h{u) + + f (1 - u 2 ) . (62) 

Taking the limit u — >• 0, we solve for C in terms of H® and H® using (|59j) : 

C = ; — — . (63) 

iu 5 - |g| 



To get the two-point functions, we need to isolate the terms in the action proportional to 

S = ^iV 3 T 6 J du d«x± [H> 2 - fH' 2 + . . .] . (64) 



H' 2 and H' 2 : 

>6^3 
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A subtle point can be made about the Gibbons-Hawking term here. In order to isolate these 
H' 2 terms in the action, we have integrated by parts terms of the form H a H'^. During this 
integration by parts, boundary terms of the form H a H' a appear at u = 0. One might think 
that these boundary terms will affect the overall normalization of the two-point functions, but 
they can't. The reason is that the Gibbons-Hawking term is precisely of a form that cancels 
these particular boundary contributions. Recall that the extrinsic curvature is defined to be 

K = -V% (65) 

where n M is a unit vector, normal to the boundary u = 0. There is a piece of K that looks 
like 



2^g^K = -V=99 uu 9 aP 9ap,u\ u=0 + • • • , (66) 
which precisely cancels the boundary contribution from terms of the form H a H'^. 

From this quadratic piece of the action, we read off the retarded Green's functions: 

2 5 vr 2 iV 3 Ty 

&ty ' ty ~ 3«(iu; - Dq 2 ) ' (67) 
2W>r B o; g 

&ty > xy ~ -&(hj-Dq*) ' (68) 
= 2Wg^ 

xy ' xy &(iu - Dq 2 ) ■ 1 ' 



where the diffusion coefficient is 



D = -A; • (70) 
4vrT v ; 



There exists a Kubo formula (see for example [7j) for these two-point functions that will 
let us calculate the viscosity rj. In particular 



7] = — lim 



q ^o q 2 y ' y 3 6 



2 5 7T 2 

- N 3 T 5 (n) 



There is a relation between the shear viscosity ()71|) and the diffusion constant (|70|): 
D = rj/(e + P), where e is the energy density and P is the pressure. The relation is a 
consequence of the conservation condition on the stress-energy tensor, d a T a ^ = 0, along 
with a linearized hydrodynamic equation for the purely spatial parts of the stress-energy 
tensor 1 

1 The bulk viscosity in these M-brane models is zero. 
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where d is the dimension of AdSd- Putting these two equations together, one finds a diffusion 
equation for the shear modes, i.e. the fluctuations of the momentum density T° l such that 
diT 0i = 0. The diffusion equation for T ty is d t T ty = Dd 2 x T ty where D = rj/(e + P). 

This relation between 77 and D will let us check that our calculations are internally 
consistent. From the thermodynamic relation between the entropy density and the pressure, 
S = dP/dT, one can calculate P. Because the stress energy tensor is traceless, for the 
M5-branes e = 5P. One finds that indeed D = rj/6P = \/4nfT. 

4.2 Vector Modes and a Diffusion Constant for M2-branes 

The same story can be repeated almost verbatim for the M2-brane case. We take the 
analogous definition for H x and H t , being careful to raise the indices of the h ty and h xy metric 
perturbations with the appropriate non-extremal M2-brane metric instead. The linearized 
Einstein's equations for H x and H t are 

co 2 H' t + q 2 H'J = 0, (73) 
H' t '-^-^(co 2 q 2 H x + q 2 2 H t ) = 0, (74) 



HZ-^H' x + -L(u 2 q 2 H t + u 2 2 H x ) = 0. (75) 
uf Af 2 



Some small changes in the wave-vector dependent pieces of the differential equations are 
apparent, but otherwise the system is virtually identical to that of (JHHJ)- dHOJ) - We combine 
the first two equations (|75|l and (|73)l to get a single equation for H' t : 

HI" + + (wl - fql)H' t + (3 - 2f)H' t = . (76) 

We make the ansatz H' t = C(l — u)^ tU)2 ^F{u) and solve for the regular function F(u) 
perturbatively in u 2 and q 2 . 

F(u) =u 2 + iu 2 (±(u 2 - u) + ±- 2 u 2 h{u) - ^/ a (u) j + f 2 {u - u 2 ) . (77) 

Taking the limit u — > 0, we solve for C in terms of H x and H? using (I74B: 



u; 2 q 2 H° x + q 2 H? 
2 (iu 2 - \ql) 



C = T , f • 
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To get the two-point functions, we need to isolate the terms in the action proportional to 

„ 2 5 /V 



H' 2 and H' 2 : 



-iV 3 / 2 T 3 J du d 3 x^ [H' t 2 -fH' x 2 + ...} 



3 4 

From here, we read off the retarded Green's functions: 



(79) 



&ty ' ty - V^-Dq 2 ) ' (80) 
2*IW 2 T 2 uq 

ty ' xy &{iu-Dq 2 ) ' 1 } 

2^ 2 nN^ 2 T 2 uj 2 

xy ' xy " 3 3 (^-ZV) • 1 } 



where the diffusion coefficient is again 



D = ——= . (83) 
4ttT V ; 



Invoking the Kubo formula, one finds from the normalization of the Green's functions 
that the shear viscosity is 

9 3 / 2 7T 

77 = ^jV 3 / 2 T 2 . (84) 
o 

We can check that the diffusion constant D = 1/4tcT is consistent with this formula for the 
viscosity, D — rj/(e + P). Tracelessness of the stress energy tensor now implies that e = 2P. 
Using the fact that S = dP/dT, one finds that indeed D = rj/3P = 1/4ttT. 



4.3 Tensor Perturbations, M5-branes, and a Kubo Formula 

We consider a perturbation of the non-extremal M5-brane metric of the form h yz ^ with 
all other h^ u = 0. We decompose such a perturbation into Fourier modes, defining 

The linearized Einstein equation for h yz becomes 

0" - + 4iK 2 - fql)<P = o • (86) 

This differential equation can be rewritten more simply as Uh z y = 0. Thus this particular 
metric perturbation can be thought of as a massless scalar. 
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We solve this differential equation in the same way we have solved the others in this 
paper: we define = C(l — u)~ tul ^ 3 F(u) and solve for F(u) perturbatively in u 5 and q 5 . The 
result is 

7/1)1- _ / 1 \ 

(87) 



F{u) = 1 - Z -fh(u) + q\ ( -\h{u) - y^f 2 (u) 



The constant C can be re-expressed in terms of the boundary value of <fi, C = (f> /F(0). 
To calculate the two point function, we isolate the term in the action proportional to <j)' 2 : 



S 



From this expression, we find that the retarded Green's function is 

2 5 7r 2 T 5 iV 3 ( . 3 



G 



yz,yz 



3 6 



iw H t^Q" 

8vrT y 



(89) 



According to [3], there is a Kubo formula which relates this Green's function to the shear 
viscosity 77. The formula states that 



77 = — lim 



lim —lm(j 

q^O (jj 



yz,yz 



2V 
3 6 



3rp5 



■N d T 



(90) 



which matches (\71\i / 



5 A^-puzzles 

The viscosities and diffusion constants calculated here present new iV-puzzles, new questions 

about why particular N dependences arise in these M-brane theories. From the gravitational 

point of view, the N dependence comes from the gravitational coupling constant /t n . What 

is missing is a direct understanding of how a field theoretic description of a stack of iV M- 

branes gives rise to the same N dependences. While for D3-branes, the iV dependence could 

be understood perturbatively from SU(N) gauge theory, for the M-branes, no equivalent 

of SU(N) gauge theory exists (as yet) as an alternate description. To review, the diffusion 

constants for the M-branes are all N independent (pTH| EH IHH|) . The viscosities (|7T| EH) 

on the other hand have the same N dependence as the entropy (0 fTT|) of the corresponding 

2 There are still other ways of calculating the viscosity through AdS/CFT correspondence. For example, 
the authors of ^2] related the graviton absorption cross section by non-extremal D3-branes to the viscosity 
of finite temperature Af = 4 SU(N) super Yang-Mills theory. Presumably a similar calculation could be 
done for the non-extremal M2- and M5-brane supergravity solutions. 
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M2- or M5-brane theory. Note that the equation D — 77/(6 + P) relates the N dependence 
of the viscosity to that of the entropy and of the diffusion constant D. But there is still at 
least one new unexplained N dependence here. 



S v D D R 



D3-branes 


^N 2 T 3 


n N 2 T 3 


1 

4ttT 


1 

2-kT 


M5-branes 


27 f N 3 T 5 


25 3 f N 3 T 5 


1 

4ttT 


3 

8nT 


M2-branes 


A T 3/2rp2 

27 


2V2tt jy3/2y2 


1 

4ttT 


3 
4ttT 



The first M-brane iV-puzzle was the entropy puzzle, the iV 3 dependence of the M5-brane 
entropy (and also the N 3 / 2 dependence of the M2-brane entropy) [TT]. For the D3-branes, 
the N 2 dependence was more or less clearly the N 2 degrees of freedom of a SU(N) gauge 
theory ^3], but what could provide N 3 degrees of freedom? 

A less well known iV-puzzle comes from the N dependence of the normalized three-point 
functions for M-branes. For D3-branes, the three-point functions scale as iV -1 , as can be 
easily understood from 't Hooft counting. However, from we know that the three-point 
functions for M5-branes scale as N~ 3 ^ 2 while for M2-branes, as iV~ 3 / 4 . 

The original motivation for writing this paper arose precisely out of these iV-puzzles, 
a feeling that these N dependences might provide a key to better understanding M-brane 
theories. In particular, the more N dependences we know, perhaps the better our chances 
of finding some pattern and unlocking the mysteries of M-theory. 3 
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